Dynamical derealization of Majorana edge states by sweeping across a quantum critical point 
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We study the adiabatic dynamics of Majorana fermions across a quantum phase transition. We show that the 
Kibble-Zurek scaling, which describes the density of bulk defects produced during the critical point crossing, is 
not valid for edge Majorana fermions. Therefore, the dynamics governing an edge state quench is non-universal 
and depends on the topological features of the system. Besides, we show that the localization of Majorana 
fermions is a necessary ingredient to guaranty robustness against defect production. 
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I. INTRODUCTION 

When certain control parameter is varied, quantum me- 
chanical fluctuations may drive a critical change of the system 
ground state [1]. Even though such phase transitions occur at 
zero temperature, they have a profound influence on phenom- 
ena as diverse as high-T c superconductivity [2], magnetism, 
or quantum Hall effects. The potential in both fundamental 
and applied research has stimulated an outgrowing interest in 
the physical community to study the effects of the dynamical 
crossing of quantum critical points [3-30]. At such points, 
the correlation length of the system diverges, and the charac- 
teristic energy gap between the ground state and the lowest 
lying excitations vanishes. Accordingly, adiabatic evolution 
is precluded, and any dynamical quench across the critical 
point is accompanied by a production of excitations in the sys- 
tem. As a result, the final state will only be partially ordered, 
displaying a non-vanishing density of defects imported from 
the quantum disordered into the ordered phases. Such density 
of defects is uniquely determined by the universality class of 
the system, and can be accurately described by the so-called 
Kibble-Zurek scaling (KZ) [3, 4]. 

The KZ scaling was formerly proposed as the mechanism 
underlying topological defect production in a cosmological 
scenario [3-5], or in classical phase transitions occurring at 
finite temperature[6, 7]. The extension of the KZ mecha- 
nism to the quantum domain was proposed in [8], and nu- 
merically confirmed in the transverse Ising model, a corner- 
stone in the theory of quantum phase transitions. This re- 
sult paved the way to a considerable amount of works dealing 
with defect production in the zero-temperature regime [9-17], 
which showed that the KZ scenario also holds for quantum 
phase transitions (see [18] for an experiment). Modifications 
to the characteristic scaling laws have been studied in order to 
account for non-linear quenches [19, 20], crossing of multi- 
critical points [21] (see [22] for a recent proposal on the im- 
plementation of tri-critical points in trapped ions), quenches 
across a critical region [23, 24], or along a gapless line [25]. 
The effect of classical noise on the driving parameter has also 
been studied in [26]. In [27-29], it was shown that the den- 
sity of defects scales with the rate of the quench even at finite 
temperatures in a regime of incipient criticality. 

Recently, we evidenced that the ground state topology of 



the system modifies the quench dynamics in a subtle man- 
ner [30]. In particular, we showed how single fermions bound 
to the edges of the Creutz ladder [31] induce an anomalous 
defect production as the magnetic flux is swept across a quan- 
tum critical point. In this paper, we reinforce the scenario, by 
providing a further evidence of anomalous defect production 
in the adiabatic dynamics of edge states. We consider a set 
of Majorana fermions (i.e. real fermions) arranged in a one- 
dimensional chain [32, 33]. We shall describe the deviations 
of the standard KZ scaling, and identify the necessary ingre- 
dients to obtain robustness of edge states even in a critical 
region by resorting to the experience gained with the Creutz 
ladder [30]. The sensitivity to boundary conditions generi- 
cally indicates that a non trivial topological order is encoded 
into the system. This is specified by the presence of states 
localized at the boundaries of the system, the so-called edge 
states [34]. Note also that topological edge states arise nat- 
urally in a wide variety of systems, such as one-dimensional 
spin models [35], the integer and fractional quantum Hall ef- 
fects (QHE) [36, 37], and have been experimentally realized 
on topological insulators [38-40]. 

The subtle difference between trivial and topological in- 
sulators (TI) is rooted on the non-trivial topological features 
of the bulk bands of the latter, such as non-vanishing Chern 
numbers in the integer QHE [41]. TI represent an intriguing 
state of matter, where edge-transport exists even in the pres- 
ence of an insulating bulk gap [42]. These edge states are 
usually localized in the interface that separates two topologi- 
cally distinct insulators, the simplest example being the inter- 
face between an IQHE sample and the vacuum [36]. Besides, 
they connect the conduction and valence bands and allow for 
conduction along the interface. Other interesting TI are the 
anomalous half-integer QHE in the honeycomb [43] or square 
lattices [44], and those ones characterized by Z2 topological 
invariants [45-50] which lead to the quantum spin Hall ef- 
fect. Depending on the symmetries of the Hamiltonian and 
the dimension of the system, TI can be classified in a periodic 
table [51, 52]. Accordingly, one of the simplest TI is the one- 
dimensional chain of Majorana fermions [32], where edge 
states correspond to bound Majorana fermions localized at the 
chain ends. The possibility to isolate Majorana fermions [53] 
represents an important step towards topological quantum 
computation [54, 55], and can in principle be achieved in the 
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vortex core of two-dimensional p + ip superconductors [56]. 
A recent proposal is to create, fuse and transport Majoranas 
by the proximity effects in tri-junctions of superconductor- 
topological insulator-superconductor [57, 58]. Therefore, the 
results on the quench dynamics of the idealized Majorana 
chain discussed in this work might have applications in the 
more realistic scenario described in [57-63]. 

This paper is organized as follows. In Sec. II we describe 
the nature of Majorana edge states in a one-dimensional lat- 
tice system with broken U(l) gauge symmetry. In Sec. Ill, 
we present the results for the quench dynamics of the Majo- 
rana chain across a quantum phase transition. In particular, 
we show that bulk states fulfill the Kibble-Zurek prediction, 
whereas topological edge state provide a clear anomaly. Fi- 
nally, we present the conclusions of this work in the light of 
topological band insulators in Sec. IV. In Appendix A, we de- 
scribe the main properties of bulk states in the periodic chain. 

II. THE MAJORANA CHAIN 

In this section, we review the properties of the Majorana 
chain [32], a fermionic one-dimensional system that breaks 
time-reversal and U(l) gauge symmetries, and presents Ma- 
jorana fermions bound to its edges. We present a detailed de- 
scription of the different phases where edge states occur, and 
discuss the nature of the quantum phase transition connecting 
the aforementioned topological phases. A discussion on the 
bulk spectrum is presented in Appendix A. 

Let us consider a system of spinless fermions hopping in 
a one-dimensional chain according to the following Hamilto- 
nian 

H = Y^(^-waja j+l +Aaja j+l -^a]aj + h.cj , (1) 

where aj (at) represent spinless fermionic annihilation (cre- 
ation) operators satisfying canonical anticommutation rela- 
tions {aj,a\} = ajal + a|a ; - = 8jk- As a second quantized 
Hamiltonian, Eq. (1) describes a mean-field Bardeen-Cooper- 
Schrieffer (BCS) superconductor [64], where w is the hopping 
amplitude, A = |A|e' e stands for the superconducting gap, and 
ji stands for the chemical potential. Due to the U( 1 ) symmetry 
breaking inherent to superconductors (i.e. absence of fermion 
number conservation), it is more appropriate to describe the 
system in terms of Majorana fermions 



which are hermitian operators ci^j-i = C2/-i)c|j = cy sat- 
isfying the Majorana fermionic algebra {e/,Cfc} = djk- Al- 
though Majorana excitations are commonly paired to consti- 
tute a standard fermion, very recent research activity proposes 
to isolate and manipulate them, being topological quantum 
computation the major motivation for doing so [54, 55]. Here, 
V = | fractional quantum Hall effect [65-67] or time-reversal 
topological insulators [57] are the most promising candidates. 



In order to understand how the Majorana chain supports 
unpaired Majorana fermions at the boundaries, we need to 
rewrite the Hamiltonian of Eq. (1) in the Majorana picture 

H = l£(tC2jC2j+l+UC 2 j-lC2j+2 + VC2j-lC2j) , (3) 

where t = w + |A|, u = — w+ |A|, and v = — ji represent the 
different couplings between Majoranas, which might be inter- 
preted as simple links in the ladder configuration of fig. 1 . 
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FIG. 1 : Scheme of the Majorana chain represented as a virtual ladder. 
The initial fermionic chain is mapped onto the Majorana fermions by 
means of Eq. (2), which then can be represented in a virtual ladder, 
where the ;-th rung contains two Majoranas c 2 j-\,c 2 j and corre- 
sponds to the y'-th site of the original chain. In the ladder scheme, the 
Majorana couplings in Eq. (3) become simple vertical and diagonal 
links. 

We shall be interested in the following regimes: 

a) Left-handed regime: This regime occurs at (t,u,v) = 
(0,2|A|,0), and presents the following zero-energy edge 
states, completely localized within the boundary rungs 

ii>/ = -oo) - C2\a), \t), = \n- ■::) = c^m, w 

where |£2) stands for the Bogoliubov vacuum (see Eq. (A7) 
in Appendix A). As represented in fig. 2(a), this is the only 
possibility to obtain localized Majoranas, since the remaining 
bulk eigenstates consist of paired of Majoranas with energies 
Ej± = ±2|A|, and eigenstates 

\£j±), = loo • • -o • • -o o> = ^(±kty-l + Cy +2 )\a). (5) 

The bulk solution consist of a flat band of plaquettes (i.e. two 
rungs) that contain a couple of Majoranas dimerized into a 
standard complex fermion. 

b) Right-handed regime: In this regime (t,u,v) = 
(2|A|,0,0), we obtain the following zero-energy edge states 
completely localized within the boundary rungs 

\i) r =\:z---° ° )=c 1 \£i), \i)r=\° i-i:)=c 3L \a). (6) 

In fig. 2(b), we represent the localized Majoranas, and the bulk 
plaquette-like eigenstates with energies E/± = ±2|A| become 

\£j±)r = |oo-"o«o-"oo) = 75(±iC2j + C 2 ;+l)|ii>. (7) 

Note how the bulk dimers containing a couple of paired Ma- 
joranas have been tilted with respect to the left-handed regime 
in fig. 2(a) 
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c) Topologically-trivial regime: In the regime \t\, \u\ <C |v|, 
every Majorana is paired up with its rung partner, and form 
plaquette-like eigenstates in flat bands £/± = ±jtl 

\e j ±) = \Z---l\l---ll)j i {Tic2 ] -i + C2 j )\eL). (8) 

Therefore, in this trivial regime edge states are forbidden (see 
fig. 2(c)). 




(a) Left-handed regime 




(b) Right-handed regime 



Iliilil 

(c) Topologically-trivial regime 

FIG. 2: Scheme of the unpaired Majorana fermions bound at the 
edges, and the dimerized Majoranas at the bulk, (a) In the left-handed 
regime, only the left-tilted diagonal links survive, and thus a couple 
of unpaired Majoranas appear at sites 2,2L— 1 of the ladder, (b) In 
the right-handed regime, only the right-tilted diagonal links survive, 
and thus a couple of unpaired Majoranas appear at sites 1,2L. (c) 
In the topologically-trivial regime, every Majorana fermion is paired 
within a single rung, and thus no edge can exist. 

We have thus shown that the Majorana chain Hamilto- 
nian introduced in Eq. (2) entails different phases where edge 
states are pinned to the system boundary, but also trivial 
phases where unpaired Majoranas cannot exist. We comment 
that the Hamiltonian in Eq. (1) could be realized with junc- 
tions involving s-wave superconductors and topological insu- 
lators [57]. In particular, these different phases can be engi- 
neered in superconductor-topological insulator tri-junctions, 
where the key parameter that determines the existence of un- 
paired Majoranas is the relative phase of the superconductors. 
In the following section, we describe the quantum phase tran- 
sitions between the phases introduced above, and study how 
the edge and bulk states behave when the system is driven 
across a quantum critical point. 



III. QUENCH BETWEEN EDGE STATES 

In this section, we study the quantum phase transition be- 
tween the left- and right-handed phases (see a) and b) above). 
Such a quantum phase transition can be characterized by t = 
and a relative parameter £ = w/|A|, so that the left-handed 
regime is realized at £ = — 1 and the right-handed at £ = 1 . In 
fig. 3, we depict the energy spectrum of the Majorana Hamil- 
tonian in Eq. (3) for £ G [—1,1], which clearly shows a quan- 
tum critical point at £ c = where the energy gap vanishes. 
Note also how two zero-energy modes, corresponding to the 



unpaired Majoranas introduced above, are present for arbi- 
trary £. Let us remark here that such edge modes are absent 
in the periodic chain (see Eq. (A6) and fig. 7 in Appendix A). 




-1 -0.5 0.5 1 

C 

FIG. 3: Energy spectrum b/\A\ of the open Majorana chain for t = 0, 
and relative parameter £ = w/\A\ £ [—1, 1]. 

The critical properties of this quantum phase transition are 
characterized by the critical exponents V = z = 1, which de- 
scribe how the coherence length |; diverges and the energy 
gap Ae vanishes at criticality <fj; ~ Ag~' ~ |£ — £ c | -1 . We 
shall consider the adiabatic dynamics of a state initially local- 
ized within the bulk or edge of the system for w(to) = — |A|. 
At the end of the quench vv(ff) = — |A|, the system shall be 
found in a state l'F(ff)) given by the superposition of the 
adiabatically-connected state, which is still localized, and de- 
localized defects which have been excited close to the critical 
point. Since we are concerned with single-particle edge- or 
bulk-states, we shall refer to the probability to create a de- 
fect Pdef = Le>o I (E^itf) | 2 , with (E\ being the eigenvectors 
of the Hamiltonian. Interestingly enough, for a linear adia- 
batic quench = — 1 + vgf with rate vq <C 1 and lasting for 
t € [0,2/vq], the Kibble-Zurek mechanism predicts that such 
probability scales as 

dv 

~ v|p , (9) 

where d is the dimension of the system [3, 4]. In the present 
phase transition, the KZ mechanism predicts a scaling ~ 
y/vQ- In the following sections, we shall confront the KZ pre- 
diction with the exact dynamics of the Majorana chain. 

A. Quench dynamics in the periodic chain 

In Appendix A, we describe the energy spectrum of the 
translationally-invariant Majorana chain, and show how the 
single-particle energy levels can be expressed as excitations 
over the Bogoliubov vacuum (i.e. yj ± |£2)). Here, we study 
the quench dynamics of the periodic chain when the parame- 
ter £ is adiabatically modified across the corresponding criti- 
cal point. In the translationally-invariant case, no edge states 
exist, and the only localized state in the negative energy band 
E = — 2|A| at £ — —1, is a plaquette-like fermion 

1^(0)) = ^(-kay_i+C2y + 2)|fl) =JV P &(f0)|a>, ( 10 > 

where c q (t ) = ^(e~T (-1 + e 1 '"), -e+T (1 + e 1 '")) deter- 
mines such initial state in momentum Nambu-representation. 



4 



Remarkably, the whole quench dynamics can be expressed 
as a collection of uncoupled two-level systems, each associ- 
ated to a Nambu spinor l P 9 . Therefore, the dynamics under 
£(/) = — 1 +VQt entails an ensemble of two-level Landau- 
Zener processes [11, 69, 70], and the final excitation proba- 
bility (compare to Eq. (12)) can be simply calculated as 



„C =1 

'q+ 



? ('/)| 2 



(ID 



where c q (tf) is the final state after the quench evolution, and 

c*~^ correspond to the positive-energy eigenstates at £ = +1. 
The evaluation of the excitation probabilities is simpler in 
this translationally invariant case, and we can thus treat larger 
chains and slower quenches. In fig. 4, we observe that the 
scaling fully agrees with the KZ-prediction P^ ~ y/vQ- 
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FIG. 4: Scaling of the defects produced by adiabatically crossing 
the Q = critical point for the periodic lattice with number sites 
L e {40, 80, 160, 320, 640, 1280, 2560}. 

In this section we have considered a single-particle initial 
state localized within a rung of the Majorana ladder. Let us 
note that the KZ-scaling also holds for an initial ground state, 
which consists of all the negative-energy modes being occu- 
pied. For such state, one usually refers to the density of de- 
fects «£j ~ which is finite in the thermodynamical limit. 
However, since we are interested in the comparison of initial 
single-particle bulk-, or edge-states, we shall refer to the prob- 
ability to create defects Pd e f rather than to a thermodynamical 
density. 



B. Quench dynamics in the open chain 

The objective of this section is to test the prediction for an 
initial plaquette-like state in the bulk, and for an initial edge 
Majorana fermion (see figs. 5(a) and 5(b) for the adiabatic 
evolution of such initial states). 

Adiabatic quench for an initial plaquette: In this case, the 
initial state corresponds to a negative-energy e = — 2|A| local- 
ized in the bulk |»P(0)) = -^(-icL-i +c L+2 )|n) for £ = -1. 
We evaluate the defects produced close to the critical point as 
follows 



Pdef=£K«f)>| 2 

£>0 



Let us emphasize that the initial bulk fermion can be ex- 
pressed as linear combination of every periodic eigenstate in 
the negative-energy band (see Eq. (A3)). In particular, the 
modes responsible for the universal KZ-scaling (i.e. those 
around the gapless mode q ~ 0) shall be populated, and one 
thus expects the KZ prediction to hold. Indeed, the numerical 
results in fig. 6(a) are in clear agreement with the predicted 
scaling P^ ~ ^JVq at the thermodynamical limit L — > °° (see 
also fig. 4 for the KZ-scaling in a periodic chain). 

Adiabatic quench for an initial edge Majorana: In this 
case, the initial state corresponds to a zero-energy Majorana 
fermion bound to the left edge ^(O)) = |1); = c 2 |£2) for 
£ = — 1. Let us note that this state also crosses the critical 
point at the gapless mode energy, and thus one would ex- 
pect a similar KZ-scaling to hold P^ ~ ^Jvq. Note how- 
ever, that the obtained scaling (fig. 6(b)) completely disagrees 
^def ~ ( v 2)°> an( l offers thus further evidence of the anoma- 
lous defect production for localized edge states [30]. The de- 
parture from the universal KZ-scaling is drastic for such Ma- 
jorana fermions. Besides, it predicts that the edge state robust- 
ness is lost as the critical point is crossed (i.e. Pdef = \ both 
for positive- and negative-energy bands). The underlying rea- 
son for such a departure is that the edge state at the critical 
point can be expressed as a linear combination of the gapless 
modes q = {0, 7f}, and thus fuses into the bulk bands. Indeed, 
we show in the next section that in the thermodynamical limit 



|l(Cc)>, ~ (c 2 + c 6 + • • • + c 2L _ o \£l) ~ -j= + yV) \Q) 

(13) 

where yj ± are the fermionic creation operators of positive- 
and negative-energy solutions associated to the mode q (see 
Eq. (16) ). Therefore, the initially unpaired Majorana at the 
left edge, becomes completely delocalized along the upper 
chain of the rung (see fig. 5(b)) at criticality £ e = 0. Being 
uniformly delocalized, it overlaps with bulk excitations with 
well-defined momenta q = {0, n}, and therefore completely 
fuses in the continuum bands. In this regard, the Majorana 
character of the initial state is completely lost across the quan- 
tum critical point, and one cannot further transport the Majo- 
rana to the right-hand side (i.e. |r) r ). 



o o . ^fc . . ^ . " & '. ' 



(a) Adiabatically connected bulk plaquettes 



o o o 



» M J> 



(12) 



(b) Adiabatically connected edge Majoranas 

FIG. 5: Adiabatically connected states for £ € [—1,1] (a) In the chain 
bulk, left- and right-tilted plaquettes with energies e = — 2|A| are 
connected |£/_); — ► \Sj-) r . (b) In the chain edges, unpaired Ma- 
jorana fermions at different edges are connected |1); —¥ |r) r , and 
|r), -> |1>, [68]. 
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Due to the finite number of modes coupled to the edge 
states, this anomalous dynamics might be viewed as a satu- 
ration effect, such as that one studied in the XY model [26]. 
Nevertheless, we remark that the effect presented in this work 
traces back indeed to the topological order of the system, since 
it only occurs for the unpaired Majorana edge states. In a triv- 
ial insulator, where in-gap edge states are simply absent, we 
should have restricted to bulk plaquette-like excitations which 
do not lead to any anomaly in KZ scaling. 



on the Creutz ladder highlights the resilience of edge states 
due to a quantum interference phenomenon that allows them 
to remain well-localized at the boundaries of the lattice. Con- 
versely, Majorana edge states become completely delocalized 
at the critical point, and this leads to a higher sensitivity at 
criticality P^f ~ Vg. Therefore, we can draw the conclusion 
that edge state robustness is non-universal, and is intimately 
related to edge state localization. 
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(a) Produced defects for a bulk state quench 
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(b) Produced defects for an edge state quench 

FIG. 6: Scaling of the defects produced by adiabatically cross- 
ing the £ c = critical point for lattice with number sites L £ 
{10,20,30,40,50,60,70,80,160}. (a) Defects P def for an initial 
negative-energy plaquette localized at the bulk, (b) Defects P^f for 
an initial zero-energy Majorana localized at the left edge of the sam- 
ple. 

Let us now compare these results with the edge-state defect 
production in the Creutz model, a ladder of usual fermions 
pierced by a perpendicular magnetic field [30]. In this sys- 
tem, a quench of the magnetic flux drives the fermions across 
a quantum phase transition of the same universality class 
v = z = d = 1 . Similarly to the Majorana chain, the defect 
production for initial bulk states fulfills f^ef ~ y/vQ in accor- 
dance to the KZ-scaling. This result emphasizes the fact that 
both models, being microscopically different, belong to the 
same universality class. Conversely, edge states in the Creutz 
ladder fulfill P^ef ~ v g 35 ' which completely deviates from the 
KZ prediction, and from the defect production of a Majo- 
rana chain P^f ~ Vg. This result gives further evidence of the 
conjecture on the non-universality of the edge state dynamics 
across a critical point, as raised in [30]. Moreover, the result 



C. Adiabatic path for edge states 

Here, we show how the unpaired Majoranas, initially local- 
ized at the chain edges, become completely delocalized over 
the whole chain when the critical point £ r = is achieved. 
Besides, we show that such delocalized states overlap with 
the continuum modes in the thermodynamic limit L — > °°, and 
thus justify the scaling P^ef ~ j( v q)°- Remarkably, the zero- 
energy modes associated to unpaired Majoranas can be ana- 
lytically obtained in the whole regime — 1 < C < 0, r/ = 0, 
and read as follows 



IKC)> = ^(« 



L+l 

■/■Teat £2 



KC)> = ~jy \C2L-l+rC2L-5- 



(14) 



where r = (1 + £)/(l — Q, and Jf^ tf are normalization con- 
stants. Note that the exponential tail of the edge states fulfill 
< r < 1 in the regime of interest. In particular, at the ini- 
tial point £ = — 1, one can easily check that these edge states 
become completely localized at the edge rungs and coincide 
with Eqs. (4). Conversely, at criticality £ c = 0, the Majoranas 
become uniformly delocalized along each chain forming the 
virtual rung 



■C6- 



(15) 



|r(«) 



( Cl +c 5 + ---+c 2L _,)|n). 



Its precisely at this critical point where they overlap with the 
bulk gapless modes (see Eq. (A3)), and thus 



rl 



HQ) 



(16) 



in) 



Consequently, the edge states will not follow the adiabatic 
path that would connect them to the edge Majoranas at £ = 
+ 1, but rather fuse into the positive- and negative-energy bulk 
bands. This explains the anomalous scaling P^ef = \ f° r edge 
states (note also that the probability to decay into the negative 
energy band would be P neg = i). 



IV. CONCLUSIONS 

In this article, we have described the quench dynamics of 
a ID-lattice model with unpaired Majorana fermions bound 
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to the boundaries [32]. The transport of Majorana fermions 
between the edges of the chain is hampered by the loss of adi- 
abaticity at the quantum critical point. The main conclusion 
drawn from this work is two-fold: edge-state dynamics across 
a critical point is not universal, and the robustness/fragility of 
the latter relies on their localization at criticality. 

We have found that the scaling of defects deviates from the 
usual Kibble-Zurek paradigm, and is non-universal. Indeed, 
the Creutz ladder [30] and the Majorana chain present dif- 
ferent scalings within the same bulk universality class. With 
respect to our results on the Creutz ladder [30], we conclude 
that such deviations are due to the localization dynamics of 
edge states. We note, in particular, that Majorana edge states 
produce a peculiar scaling of excitations due to their interest- 
ing derealization dynamics. At the quantum critical point, 
Majorana edge states 'fuse' with the bulk excitations, and 
thus become delocalized over the system. This should be 
contrasted with the phenomenology evidenced for the Creutz 
ladder, where edge states stay localized for any value of the 
control parameter. Such a difference, in turn, traces back 
to the physical content of such states: for the Creutz ladder, 
standard fermions are pinned to the ladder edges due to in- 
terference effects caused by additional magnetic fields; for 
the Majorana chain, they are induced by the U(l) symmetry 
breaking term provided by a superconducting gap. Despite of 
this quantitative difference in the microscopic origin, we note 
that both systems belong to the class of topological band in- 
sulators [51, 52]. Therefore, we may conclude that quench 
dynamics of boundary states in topological insulators should 
provide additional anomalies in the defect production rate. 
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Appendix A: Bulk spectrum of the Majorana chain 

In this Appendix, we describe the bulk spectrum of the Ma- 
jorana chain for periodic boundary conditions (ol+i = a\). 
Let us first point out that the Hamiltonian in Eq. (1), after a 
standard Jordan- Wigner transformation [71], describes a set 
of localized spins governed by 

H = HA ' COS0 t (C sec0 + 1 W+i + (CsecG - 1)ctJctJ +1 

Z 7=1 

+ tan 6 {&ja> j+i + a]aj +i ) + h(a] + 1 ) , 

(Al) 

where h = jU/|A| cos 9 plays the role of an external transverse 
magnetic field, and the superconducting phase leads to the 
local lattice anisotropy in the exchange couplings, namely, 
J x = (£sec0 + l), and J y = (£sec0 - l), with £ = w/\A\. 
Notice also the appearance of additional spin coupling terms 
OjOj +1 , which are absent in standard anisotropic XY mod- 
els [72, 73]. Accordingly, a spin picture offers an alternative 



perspective on the U(l) symmetry -broken fermion system in 
Eq. (1). In momentum space, the fermionic operators become 
aj = -jjY* q a q zT lqj ', where q € [—n,n] lies in the Brillouin 
zone and we have assumed unit lattice spacing. The initial 
Hamiltonian in Eq. (1) can be expressed as follows 



ff = £M*ff,*P„ H q = 

q 




(A2) 



where we have introduced the Nambu spinor m q = (a q ,a_ q )', 
and £ q = — (wcosq + ju/2), A q = i|A|e~' e sing. The energy 
spectrum is obtained after a Bogoliubov transformation is 
performed, in a similar procedure as the Bardeen-Cooper- 
Schrieffer (BCS) Hamiltonian is diagonalized [64]. Trans- 
forming the fermion operators according to 

Yq+ = u* q a q + v q ai q , y q - = -v*a q + u q a_ q , (A3) 

one obtains the energies e q ± = ±e q = ±^J Q + lA^I 2 corre- 
sponding to the positive- and negative-energy solutions y q ±, 
defined through the following parameters 

1 r~T q ie- i9 sgn( g ) [~T q 

Besides, the solutions fulfill J- q ± = y qzf , which allows us to 
restrict the Hamiltonian to half-Brillouin zone 

H = E ( l£ d+ Ki+ - 2e «i- r«- )■ ( A5 > 

Therefore, the final energy spectrum for a periodic Majorana 
chain is 



E q ± = ±2e, 7 = ±2|A|^(77 + Ccos<7) 2 + sin 2 <7, ( A 6) 

where we have introduced the relative parameters £ = 
w/|A|,tj = ju/2|A| used to study the quantum phase transi- 
tions in the open chain. In fig. 7, we represent the energy 
spectrum as a function of £ for 77 = 0. We clearly identify 
a non-topological two-band insulator since no mid-gap edge 
states occur at zero-energies. This fact should be compared 
to the open ladder figs. 3, where zero-energy edge states lo- 
calized at the interface between the topological insulator and 
vacuum arise. In these figures, it is also clear that the en- 
ergy gap vanishes at the critical points £ c = 0, indicating thus 
a quantum phase transition. Furthermore, from the exact en- 
ergy dispersion in Eq. (A6), one can show that the critical 
exponents of both transitions are z = V = 1 . 

Let us finally comment on the system vacuum and single- 
particle states. One may define the Bogoliubov vacuum as 
a state satisfying y c/ ±|£2) = 0,Vg € BZ. This vacuum can be 
constructed from the Fock state vacuum |0) as follows 

i«> = n^+r,/-io> = nK+ v A4)i°)' < A7 > 

9 q 

which shows that the Bogoliubov vacuum contains the usual 
Fock vacuum, where a number of exotic spinless Cooper pairs 
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have been condensed into. One can check that in the limit of 
vanishing gap (A — > 0), the standard Fock vacuum is recov- 
ered. The positive- and negative-energy single-particle levels 
can be built from the Bogoliubov vacuum as = yLjil). 
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